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Abstract We have constructed a new type of amplifier whose primary purpose is the read- 
out of superconducting quantum bits. It is based on the transition of an RF- 
driven Josephson junction between two distinct oscillation states near a dynam- 
ical bifurcation point. The main advantages of this new amplifier are speed, 
high-sensitivity, low back-action, and the absence of on-chip dissipation. Us- 
ing pulsed microwave techniques, we demonstrate bifurcation amplification in 
nanofabricated Al junctions and verify that the performance predicted by theory 
is attained. 

1. Introduction 

Josephson first noted that the superconducting tunnel junction can be viewed 
as a non-linear, non-dissipative, electrodynamic oscillator [1]. We exploit this 
non-linearity to produce a new type of high-sensitivity amplifier, the Josephson 
Bifurcation Amplifier (JBA). No shunt resistors are required in our amplifica- 
tion scheme, and it is thus possible to take advantage of the elastic charac- 
ter of the junction and eliminate on-chip dissipation, thereby minimizing the 
back-action of the amplifier. The combination of high-sensitivity and minimal 
back-action makes the JBA well-suited for measurements on quantum systems 
such as superconducting qubits, and make it a strong candidate for reaching 
the quantum noise limit. 

The operation of the JBA is represented schematically in Fig. 1 . The cen- 
tral element is a Josephson junction, shunted with a lithographic capacitor, 
whose critical current In is modulated by an input signal (input port). Cou- 
pling between the junction and the input signal can be achieved through dif- 
ferent schemes, examples of which involve placing the JBA in a SQUID loop 
[2] or in parallel with a SSET [3]. The junction is driven with a pure AC 
signal irf sin(cjt) in the microwave frequency range fed via a transmission 
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line through a circulator (drive port). In the underdamped regime, for certain 
values of oj and irf, two possible oscillation states which differ in amplitude 
and phase (denoted "0" and "1") can coexist. The reflected component of the 
drive signal, measured through another transmission line connected to the cir- 
culator (output port), is a convenient signature of the junction oscillation state. 
At the bifurcation point where switching between oscillation states occurs, the 
system becomes infinitely sensitive, in the absence of thermal and quantum 
fluctuations, to variations in 1$. At finite temperature, the energy stored in the 
oscillation can always be made larger than thermal fluctuations by increasing 
the scale of Iq, thus preserving sensitivity. Small variations in Iq are trans- 
formed into readily discernible changes in the escape rate Toi from state to 

1. Back-action is minimized in this arrangement since the only fluctuations 
felt at the input port arise from the fluctuations of the 50 Q drive port whose 
dissipative elements are physically separated from the junction via a transmis- 
sion line of arbitrary length and can therefore be thermalized efficiently to base 
temperature. Additionally, the frequency band over which the back-action con- 
tributes is narrow, and well controlled. 

In section 2, simplified expressions adapted from the theory of activated 
escape in a driven non-linear oscillator [4] are presented. Details of the devices 
and the measurement setup are presented in Section 3. Experimental results are 
given in Section 3, and concluding remarks are in Section 4. 

2. Theory 

The tunnelling of Cooper pairs manifests itself as a non-linear inductance 
that shunts the linear junction self-capacitance Cj, formed by the junction 
electrodes and the tunnel oxide layer. The constitutive relation of the non- 
linear inductor can be written as I(t) = IosinS(t), where I(t), 5 (t) = 
fli^dt'Vit^/ipa and V(t) are the current, gauge-invariant phase-difference 
and voltage corresponding to the inductor, respectively, while the parame- 
ter Iq is the junction critical current. Here ipo = h/2e is the reduced flux 
quantum. The dynamics of the junction are given by the time evolution of 5, 
which exhibits the motion of a phase particle in a cosine potential U(5) = 
—(PqIqcos(5). For small oscillation amplitude about the potential minima, 
the frequency of oscillation is given for zero DC bias current by the plasma 
frequency lopq = 1/y/LjCj where Lj = po/Io is the effective junction in- 
ductance. As the oscillation amplitude increases, the potential "softens" and 
oop decreases, an effect which has been measured in both the classical and 
quantum regime [5-8]. A more dramatic non-linear effect manifests itself if 
the junction is driven with an AC current irf smut at a frequency to slightly 
below LOpQ. If the quality factor Q = Cjupo/ Re[Z~ 1 (ujpo)] is greater than 
y/3/2a, where Z(upo) is the impedance of the junction electrodynamic envi- 
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ronment and a = 1—lu/ujpo the detuning parameter, then the junction switches 
from one dynamical oscillation state to another when irf is ramped above a 
critical value Ib [9]. For ipp < Ib, the oscillation state is low-amplitude 
and phase-lagging while for ipp > Ib, the oscillation state is high-amplitude 
and phase-leading. This generic non-linear phenomenon, which we refer to 
as "dynamical switching", is reminiscent of the usual "static switching" of the 
junction from the zero- voltage state to the voltage state when the DC current 
bias exceeds the critical current Iq [10]. However, an important distinction be- 
tween dynamical and static switching is that in dynamical switching, the phase 
particle remains confined to only one well of the junction cosine potential, and 
the time-average value of 5 is always zero. The junction never switches to 
the voltage state, and thus no DC voltage is generated. Also, for dynamical 
switching, the current Ib depends both on Q and on the detuning a. 

In presence of the microwave drive ipp sin(u;t), the oscillations in the junc- 
tion phase can be parameterized using in-phase and quadrature phase compo- 
nents 5(t) = Su sm(ut) + S± cos(ujt) (higher harmonics of oscillation are neg- 
ligible). The two oscillation states appear as two points in the (#||,<5j_^ plane 
and are denoted by vectors labelled and 1 (see Fig. 2). The error-current [11] 
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Figure 1. Schematic diagram of the 
Josephson bifurcation amplifier. A junc- 
tion with critical current Io, parametrically 
coupled to the input port, is driven by an 
RF signal which provides the power for 
amplification. In the vicinity of the dynam- 
ical bifurcation point inp = Ib, the re- 
flected signal phase <f> depends critically on 
the input signal. 



Figure 2. Poincare section of an 
RF-driven Josephson junction in the 
bistable regime (a = (1 — ui/uj p ) — 
0.122, i RF /lB = 0.87). The two stable 
oscillation states, labelled by and 1, are 
indicated by white line segments. Point S 
which lies on the separatrix is the saddle 
point at which the escape trajectory from 
state (dashed line) meets the retrapping 
trajectory into state 1 (solid line). 
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which describes the generalized force felt by the system is also plotted in the 
(<5j| , <5_l) plane. Its value goes to zero at the attractors corresponding to states 
and 1 and also at a third extremum which is the dynamical saddle point. Also 
shown in Fig. 2 is the calculated escape trajectory [12] from state (dashed) 
and the corresponding retrappping trajectory [13] into state 1 (solid line). Fig. 
2 has been computed for a = 0.122, Q = 20 and irf/Ib = 0.87. These 
values correspond to typical operating conditions in our experiment. The dy- 
namical switching from state to 1 is characterized by a phase shift given here 
by tan" 1 (s^ — 5®)/(5]_ — 5j_^j = — 139 deg. Using the junction phase- 
voltage relationship and the transmission line equations, we can calculate the 
steady-state magnitude and phase of the reflected microwave drive signal. The 
change in the oscillation of S results in a shift of the reflected signal phase 
A0oi = 89 deg. Since there is no source of dissipation in the junction chip, 
there should be no change in the magnitude of the reflected signal power, even 
though ^/(<5} - <5}) 2 + (5]_ - 5° ± ) 2 + 0. 

Both static and dynamical switching can be described by an Arhennius law 
in which the escape rate Toi = (w a tt/2vr) exp(— AU/k^T) is written as the 
product of an attempt frequency io att /2tt and a Boltzman factor which contains 
the potential barrier height AU and the system temperature T. For the case of a 
DC current bias, the cosine potential, near the switching point, is approximated 
as a cubic potential with height A U st = (4\/2/3) (polo (1 — iDc/Io) 3 ^ 2 where 
idc is the bias current. The attempt frequency is the plasma frequency up. In 
the absence of fluctuations, the characteristic current at which switching oc- 
curs is Iq- For the AC driven junction, the dynamical switching from oscilla- 
tion state to 1 can be cast in a similar form using the model of a particle 
in a cubic metapotential [4]. In this case, the effective barrier height is, to 
lowest order in l/(aQ), AU dyn = u dyn (l - (i RF / 'I B ) 2 ) 3/2 with u dyn = 
64?i/(18ev / 3) a (l — a ) 3 - The attempt frequency in the metapotential is 
given by to a = co a0 (1 - (i RF / 'Ib) 2 ) 1/2 with cu a0 = 4/(3^3 RC) (u p - u) 2 . 
The bifucrcation current Ib where the state ceases to exist is given by Ib = 
16/(3^3) a 3 / 2 (l- afl 2 I . 

3. Devices and Setup 

Typical junction fabrication parameters limit the plasma frequency to the 20 
-100 GHz range where techniques for addressing junction dynamics are in- 
convenient. We have chosen to shunt the junction by a capacitive admittance 
to lower the plasma frequency by more than an order of magnitude and attain 
a frequency in 1-2 GHz range (microwave L-band). In this frequency range, a 
simple on-chip electrodynamic environment with minimum parasitic elements 
can be implemented, and the hardware for precise signal generation and pro- 
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cessing is readily available. In the first step of sample fabrication, a metallic 
underlayer - either a normal metal (Au, Cu) or a superconductor (Nb) - was 
deposited on a silicon substrate to form one plate of the shunting capacitor, 
followed by the deposition of an insulating S13N4 layer. Using e-beam lithog- 
raphy and double-angle shadow mask evaporation, we subsequently fabricated 
the top capacitor plates along with a micron sized AI/AI2O3/AI tunnel junc- 
tion. The critical current of the junction was in the range Iq = 1 — 2 fiA. By 
varying both the dielectric layer thickness and the pad area, the capacitance C 
was varied between 16 and 40 pF. 

The junction + capacitor chip is placed on a microwave circuit-board and is 
wire-bonded to the end of a coplanar stripline which is soldered to a coaxial 
launcher affixed to the side wall of the copper sample box. We anchor the RF 
leak-tight sample box to the cold stage of a 3 He refrigerator with base temper- 
ature T = 280 mK. The measurement setup is schematically shown in Fig. 3. 
Microwave excitation signals are coupled to the sample via the -13 dB side port 
of a directional coupler. The reflected microwave signal passes through the di- 
rect port of the coupler, and is amplified first using a cryogenic 1.20—1.85 GHz 
HEMT amplifier with noise temperature T/v = 4 K. A DC bias current can be 
applied to the junction by way of a bias tee. We use cryogenic attenuators, iso- 
lators, and specially developed dissipative microstrip filters on the microwave 
lines in addition to copper-powder and other passive filters [6] on the DC lines 
to shield the junction from spurious electromagnetic noise. In the first set of 
experiments which probe the plasma resonance, a vector network analyzer was 
used to both source a CW microwave signal and to analyze the reflected power 
[14]. The dynamics of the transition between the two oscillation states was then 
probed using microwave pulses [15], generated by the amplitude modulation 
of a CW source with a phase-locked arbitrary waveform generator with 1 ns 
resolution. For the pulsed experiments, the reflected signal was mixed down 
to 100 MHz and digitally demodulated using a 2 GS/s digitizer to extract the 
signal phase cf). 

4. Results 

We first probed the drive current dependence of the reflected signal phase 
{irf) by applying a 4 /xs long symmetric triangular shaped pulse with a peak 
value 0.185 Jo- The demodulated reflected signal was divided into 20ns sec- 
tions, each yielding one measurement of <j) for a corresponding value of irf. 
The measurement was repeated 6 x 10 5 times to obtain a distribution of (^{irf)- 
In Fig. 4, the mode of the distribution is plotted as a function of irf/Iq- 
For irf/Iq < 0.125, the bifurcation amplifier is always in state 0, 4> is con- 
stant and assigned a value of deg. As the drive current is increased above 
irf/Iq = 0.125, thermal fluctuations are sufficiently large to cause transitions 
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Figure 3. Schematic of the measurement setup. Thick lines correspond to 50 O. coaxial trans- 
mission lines. The network analyzer is used for CW measurements. For probing the dynamics, 
the sample is switched to the pulse generator and phase detector. A lumped element model 
for the junction chip and measurement line is shown in the lower right. The two ideal current 
sources actually represent external sources. 



to the 1 state. In the region between the two dashed lines at irf/Io = 0.125 
and irf/Io = 0.160, <\> displays a bimodal distribution with peaks centered at 
and 74 deg with the latter corresponding to the amplifier in the 1 state. The 
dotted line in Fig. 4 is the average reflected signal phase ((/>). When irf/Iq is 
increased above 0.160, the system is only found in state 1. In the decreasing 
part of the irf ramp, the system does not start to switch back to state until 
Irf/Io = 0.065. The critical switching currents Ir for the — > 1 transition 
and ijj for the 1 — > transition, calculated from numerical simulations to treat 
the inductance of wire bonds, are denoted with ticks in Fig. 4, and are in good 
agreement with experiment. The hysteresis Ig < Ir is a consequence of the 
asymmetry in the escape barrier height for the two states. Thus, the — > 1 
transition at irf = Ir is nearly irreversible, allowing the bifurcation amplifier 
to latch and store its output during the integration time set by the sensitivity of 
the follower amplifier. 

To determine the sensitivity of the bifurcation amplifier, we have character- 
ized in detail the switching in the vicinity of the — ► 1 transition. We excited 
the system with two different readout pulse protocols. In the first protocol, 
the drive current was ramped from to its maximum value in 40 ns and was 
then held constant for 40 ns before returning to 0. Only the final 20 ns of the 
constant drive period were used to determine the oscillation phase with the first 
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20 ns allotted for settling of the phase. Histograms taken with a 10 MHz acqui- 
sition rate are shown in Fig. 5. In the upper panel, the two peaks corresponding 
to states and 1 can easily be resolved with a small relative overlap of 10~ 2 . 
The width of each peak is consistent with the noise temperature of our HEMT 
amplifier. In this first method, the latching property of the system has not been 
exploited. In our second protocol for the readout pulse, we again ramp for 
40 ns and allow a settling time of 20 ns, but we then reduce the drive current 
by 20% and measure the reflected signal for 300 ns. In that latter period, what- 
ever state was reached at the end of the initial 60 ns period is "latched" and 
time is spent just increasing the signal/noise ratio of the reflected phase mea- 
surement. As shown in the lower panel of Fig. 5, the two peaks are now fully 
separated, with a relative overlap of 6 x 10~ 5 allowing a determination of the 
state 1 probability with an accuracy better than 10~ 3 . This second protocol 
would be preferred only for very precise time-resolved measurements of 1$ or 
for applications where a low-noise follower amplifier is impractical. 

A third experiment was performed to study the state 1 switching probabil- 
ity Poi (*rf) for different values of the temperature T and Jo, the latter being 
varied with a magnetic field applied parallel to the junction plane. Using the 




0.00 0.05 I- 0.10 0.15 l R 0.20 



Wo 

Figure 4. Hysteretic variation of the re- 
flected signal phase with drive current 
irf/Io- Symbols denote the mode of <f>, 
with up and down triangles corresponding 
to increasing and decreasing irf = Ib, 
respectively. The dotted line is ((j>). The 
calculated bifurcation points, /§ and Ib, 
are marked on the horizontal axis. The 
and 1 phase states are reminiscent of the 
superconducting and dissipative states of 
the DC current biased junction. 
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Figure 5. Histograms of the reflected 
signal phase (f> at irf/Io = 0.145. The 
upper histogram contains 1.6 x 10 6 counts 
with a measurement time r m = 20 ns. The 
lower panel, taken with the latching tech- 
nique, has 1.5 x 10 5 counts with a mea- 
surement time r m = 300 ns. Data here has 
been taken under the same operating con- 
ditions as in Fig 4. The dashed line repre- 
sents the discrimination threshold between 
the and 1 state. 
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first readout protocol and the discrimination threshold shown in Fig. 5, we 
obtain the switching probability curves shown in Fig. 6. Defining the dis- 
crimination power d as the maximum difference between two switching prob- 
ability curves which differ in Iq we find that at T = 280 mK, d = 57% for 
AIq/Iq = 1% - the typical variation observed in a superconducting charge- 
phase qubit [16]. The switching probability curves should shift according to 
(AI B /I B )/(AI /I ) = 3/ (4a) - 1/2 + 0(l/(aQ) 2 ), which for our case 
takes the value 5.6. In Fig. 6, the curves are shifted by 6%, which agrees well 
with this prediction. For the case of the DC current biased junction, similar 
curves would shift only by 1% since the switching current is Iq itself. Com- 
parable discrimination power using DC switching has only been achieved in 
these devices at T < 60 mK. As the temperature is increased, the switching 
probability curves broaden due to increased thermal fluctuations and the dis- 
criminating power decreases: at T = 480 mK, d = 49%. 

Finally, we determined the escape rate Toi (irf , Iq , T) as a function of zrf 
by measuring the time dependence of the switching probability, using a method 
previously applied to the determination of the static switching rates to the volt- 
age state [17]. After the initial ramp (40 ns) and settling period (20 ns), the 
reflected signal phase was extracted every 20 ns for a duration of 1 [is. By 
repeating this measurement, we generated switching probability histograms 
which we analyzed as Poi (t) = 1 — exp(— Tqi ■ t). To obtain the escape 
rate at different temperatures, two different techniques were used. In the first 
method, we varied the temperature of the cryostat and used a magnetic field 
to keep the critical current constant at Jo = 1-12 /j,A. In the second method, 
Iq was kept fixed at 1.17 fiA, and a 1 — 2 GHz white noise source irradiating 
the junction was used to increase the effective temperature. In Fig. 7 we show 
the drive current dependence of the escape rate as (m(27ra; a /Toi)) 2 / 3 plotted 
versus i 2 RF for two different sample temperatures. Data in this format should 
fall on a straight line with a slope s (T) proportional to (ud yn /kBT) 2 ^ 3 . A 
trace taken at T = 500 mK is also shown in Fig 7. 

In parallel with these dynamical switching measurements, we ran static 
switching measurements to obtain an escape temperature T s e / C . Due to in- 
sufficient filtering in our RF amplifier line outside the measurement band, 
T^ c exceeded T by 60 mK. Using uf^ and s(T) we can cast the results 
of the dynamical switching measurements into a dynamical escape tempera- 
ture T^ c n = u%jZ/k B s (T) 3/2 . We plot versus T s e / C in the inset of Fig. 
7. The agreement is very good, and only deviations at the highest temperatures 
are observed. Analyzing the dynamical switching data with Tjf c in place of 
T, we extract a value of Ud yn = 10.7 K from the T = 280 mK data with 
Iq = 1.17 [iA while the calculated value keeping higher order terms in 1/aQ 
is< n c = H.0K. 
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Figure 6. Switching probability curves 
at T = 280 mK and T = 480 mK as a 
function of the drive current irf- The dis- 
crimination power d is the maximum dif- 
ference between two curves at the same 
temperature which differ by approximately 
l%in/ . 



Figure 7. Escape rate as a function 
of drive power for two different operating 
temperatures with Io = 1.12 /iA. The 
inset shows the relationship between dy- 
namic and static escape temperatures when 
varying either the sample temperature or 
the injected noise power. 



5. Conclusion 

With the JBA operating at T^ c = 340 mK, it is possible to resolve with a 
signal/noise ratio of 1 a 10 nA variation in 1$ in a total time < 80 ns, corre- 
sponding to a critical current sensitivity of S 1 ^ 2 = 3.3 x 10~ 6 A • Hz -1 / 2 . This 
value is in agreement with the prediction from the analytical theory S 1 / 2 = 

1 /2 

V^Rf/Io, ®) (ksT/ipo) ■ Trxi , where 7] ps 1.4 near the bifurcation point and 
tpo = h/2e. The advantage of the bifurcation amplifier over SQUIDs [18] re- 
sides in its extremely low back-action. Since there is no on-chip dissipation, 
the only source of back-action is the matched isolator load, which is efficiently 
thermalized at T = 280 mK. An important point is that in the JBA, only fluc- 
tuations from the load that occur in a narrow band centered about the plasma 
frequency contribute to the back-action, whereas in the SQUID noise from 
many high frequency bands is also significant. Finally, the bifurcation am- 
plifier does not suffer from quasiparticle generation associated with hysteretic 
SQUIDS [2] and DC current-biased junctions [3] which switch into the voltage 
state. Long quasiparticle recombination times at low temperatures limit the ac- 
quisition rate of these devices while the recombination process itself produces 
excess noise for adjacent circuitry [19]. 

In conclusion, the JBA is competitive with the SQUID for applications 
where low back-action is required. Its speed, suppression of on-chip dissi- 
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pation, and latching make it ideal for the readout of superconducting qubits. 
At temperatures such that Tj^ < 60 mK, the discrimination power would be 
greater than 95%, hence permitting stringent tests of Quantum Mechanics, like 
the violation of Bell's inequalities. 
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